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1 Introduction 



In this decade exciting experimental results have been obtained in heavy baryon spec- 
troscopy. During these years, the \ + and \ antitriplet states, A+, S+, S° and A+(2593), 
S+(2790), S°(2790) and the \ + and | + and sextet states, Q*, £*, S* have been observed in 
experiments [1]. Among the s-wave bottom hadrons, only E£, and Qf, have been 

discovered. Moreover, in recent years many new states have been observed by BaBar and 
BELLE collaborations, such as, X(3872), F(3930), Z(3930), X(3940), F(4008), Z+(4050), 
F(4140), X(4160), Z 2 (4250), F(4260), F(4360), Z+(4430), and F(4660) which remain 
unidentified. 

Of course, establishing these states is a remarkable progress in hadron physics. It is 
expected that LHC, the world's largest-highest-energy particle accelerator, will open new 
horizons in the discovery of the excited bottom baryon sates [2]. The experimental progress 
on heavy hadron spectroscopy stimulated intensive theoretical studies in this respect (for 
a review see [3,4] and references therein). A detailed theoretical study of experimental 
results on hadron spectroscopy and various weak and strong decays can provide us with 
useful information about the quark structure of new hadrons at the hadronic scale. 

This scale belongs to the nonperturbative sector of QCD. Therefore, for calculation of 
the form factors in weak decays and coupling constants in strong decays, some nonper- 
turbative methods are needed. Among many nonperturbative methods, QCD sum rules 
[5] is more reliable and predictive. In the present work, we calculate the strong coupling 
constants of light pseudoscalar mesons with sextet and antitriplet baryons, in light cone 
version of the QCD sum rules (LCSR) method (for a review, see [6]). Note that some of 
the strong coupling constants have already been studied in [7-9] in the same framework. 

The outline of this article is as follows. In section 2, we demonstrate how coupling 
constants of pseudoscalar mesons with heavy baryons can be calculated. In this section, 
the LCSR for the heavy baryon-pseudoscalar meson coupling constants are also derived 
using the most general form of the baryon currents. Section 3 is devoted to the numerical 
analysis and a comparison of our results with the existing predictions in the literature. 

2 Light cone QCD sum rules for the coupling con- 
stants of pseudoscalar mesons with heavy baryons 

Before presenting the detailed calculations for the strong coupling constants of pseudoscalar 
mesons with heavy baryons, we would like to make few remarks about the classification of 
heavy baryons. Heavy baryons with a single heavy quark belong to either SU(3) antisym- 
metric 3p or symmetric 6^? flavor representations. Since we consider the ground states, the 
total spin of the two light quarks must one for 6p and zero for 3 F , due to the symmetry 
property of their colors and flavors, as a result of which we can write J p = | + /| + for 6^ 
and J p = | + for 3^. Graphically, 6p and 3^ representations are given in Fig. (1), where 
a, a + 1, a + 2 determine the charges of baryons (a = —1 or 0), and the asterix (*) denote 
J p = | + states. In this work, we will consider only J p = | + states. 

After this preliminary remarks, we proceed by calculating the strong coupling constants 
of pseudoscalar mesons with heavy baryons within the LCSR. For this purpose, we start 
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by considering the following correlation function: 

= i J d 4 xe ipx (V(q) \T {rf i) {x)ff j \{S)}\ 0) , (1) 

where V(q) is the pseudoscalar-meson with momentum q, 77 is the interpolating current 
for the heavy baryons and T is the time ordering operator. Here, % — 1, j — I describes 
the sextet-sextet, i — 1, j = 2 corresponds to sextet-triplet, and i = 2, j = 2 describes 
triplet-triplet transitions. For convenience we shall denote = nJ 12 ) = and 

n (22) =n ( 3) _ The 

sum rules for the coupling constants of pseudoscalar mesons with heavy 
baryons can be obtained by calculating the correlation function (1) in two different ways, 
namely, in terms of the hadrons and in terms of quark gluon degrees of freedom, and then 
matching these two representations. 

Firstly, we calculate the correlation function (1) in terms of hadrons. Inserting complete 
sets of hadrons with the same quantum numbers in the interpolating currents and isolating 
the ground states, we obtain 

nW) = (0 1^(0)1 B 2 (p)) (B 2 (p)V{q) I Bi(p + g)> (Bi(p + g) 1^(0)1 0) + . . . (2) 

(jp, _ _|_ qY _ m 2j 

where |-B 2 (p)) and \B 1 (p + q)) are the \ states, and m 2 and m\ are their masses, respec- 
tively. The dots in Eq. (2) describe contributions of the higher states and continuum. It 
follows from Eq. (2) that in order to calculate the correlation function in terms of hadronic 
parameters, the matrix elements entering to Eq. (2) are needed. These matrix elements 
are defined in the following way: 

(0|j7 (i) |5(p)> = A i «(p), 

(B(p + q)\r ] W\0) = \ j u{p + q) , 

(B(p)V(q) I B{p + q)) = gu{p)i l5 u{p + q) , (3) 

where A, and Xj are the residues of the heavy baryons, g is the coupling constant of pseu- 
doscalar meson with heavy baryon and u is the Dirac bispinor. 

Using Eqs. (2) and (3) and performing summation over spins of the baryons, we obtain 
the following representation of the correlation function from the hadronic side: 

n (u) = i— 2\T7~ \2 21" ($75 + other structures), 

(p z — m^) lip + qy — mf\ I ) 

(4) 

where we kept the structure which leads to a more reliable result. 

In order to calculate the correlation function from QCD side, the forms of the interpo- 
lating currents for the heavy baryons are needed. The general form of the interpolating 
currents for the heavy spin | sextet and antitriplet baryons can be written as (see for 
example [10]), 

^ = -^e abc {(qfCQ b ) l5 q c 2 + p{qfC l5 Q b )q c 2 - [ (Q aT C^) 75^ + ^ (^75^) ^] } , 
V { Q nU ~ t] = ^ abC {^(lfCqi)l 5 Q c + 2/3(gfC 7 5^)g c + {qfCQ b ) l5 q c 2 + p{qfC l5 Q b )q c 2 

+ (Q^Cq^qt + P^C^qfyl} , (5) 
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where a, ,b,c are the color indices and /3 is an arbitrary parameter. It should also be noted 
that the general form of interpolating currents for light spin 1/2 baryons was introduced in 
[11] and (3 = — 1 corresponds to the Ioffe current [12]. The quark fields qi and q 2 for the 
sextet and antitriplet are presented in Table 1. 
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Table 1: The quark flavors gi and g 2 for the baryons in the sextet and the antitriplet 
representations 

As has already been noted, in order to calculate the coupling constants of pseudoscalar 
mesons with heavy baryons entering to sextet and antitriplet representation, the calculation 
of the correlation function from QCD part is needed. Before calculating it, we follow the 
approach given in [13-17] and try to find relations among invariant functions involving 
coupling constants of pseudoscalar mesons with sextet and antitriplet baryons. We will show 
that the correlation functions responsible for coupling constants of pseudoscalar mesons (P) 
with sextet-sextet (SS), sextet-antitriplet (SA) and antitriplet-antitriplet (AA) baryons 
can each be represented in terms of only one invariant function. Of course, the form of the 
invariant functions for the couplings SSP, SAP and AAP are different in the general case. 
It should be noted here that the relations presented below are all structure independent. 

We start our discussion by considering the sextet-sextet transition, concretely. Consider 
the E° — » S°7T° transition. The invariant function for this transformation can be written in 
the following form 

n s°->£<v° = ^ n (D dj b) + ^ n 'a) dj b) + ^ n (D dj b)i (6) 

where the interpolating current of n° meson is written as 

Obviously, the relations g nau = ~9ndd = 775 , fiVe* = hold for the 7r° meson. The 
invariant functions III, 11^ and II2 describe the radiation of it meson from u, d and b 
quarks of E° baryon, respectively, and they can formally be defined as: 

nS 1} (M,rf,6) = (uu 0) , 

n'^u,^) = <Jrf|s°s°| 0) , 

U^\u,d,b) = <5fe|E°E°|0> . (7) 

It follows from the definition of the interpolating current of baryon that it is symmetric 
under the exchange u ■<-)■ d, hence n i ^(u, d, 6) = n^(d, -u, 6). Using this relation, we 
immediately get from Eq. (6) that, 



n s°^£;y>_ 1 



V2 



U?(u,d,b)-U?(d,u,b)] , (8) 
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and one can easily see that in the 67/2(2)/ limit, n 1 ^ 2 " 71 " = 0. 

The invariant function responsible for the E^" — > £^7r° transition can be obtained from 
the E° — > S°7r° case by making the replacement d — » u, and using E° = — \/2E^, from 
which we get, 

AU { ^(u,d,b) = -2{uu |E+E+| 0) . (9) 

Appearance of the factor 4 on the left hand side is due to the fact that each E^~ contains 
two u quark, hence there are 4 possible ways for radiating it from the u quark. Making 
use of Eq.(8), we get 

nE+^E+Tr" = V2U?(u,u,b). (10) 

The invariant function describing E^ — > S^7r° can easily be obtained from the E° — > 
S°7T° transition by making the replacement u — >■ d and taking into account E°(ii — > d) — 
y/2T. b . Performing calculation similar to the previous case, we get 

n E ^ s ^° = V2u?(d,d,b). (ii) 

Now, let us proceed to obtain the results for the invariant function involving S 6 — > 
S 6 transition. The invariant function for this transition can be obtained from the 

E£ ->■ S^7T° case using the fact that E b ° = E£(d ->• s) and Z b ~ — E£(u ->• s). As a result, 
we obtain 

n s ^ s ^° = --^=nS 1) (rf, s ,6). (12) 

Obtaining relations among the invariant functions involving charged 7r =l= mesons requires 
more care. In this respect, we start by considering the matrix element (dd |E°S°| 0), where 
d quarks from the S° and S° from the final dd state, and u and b quarks are the spectators. 
The matrix element (ud |E^£°| 0)describes the case where d quark from S° and u quark 
from form the ud state and the remaining u and b are being again the spectators. One 
can expect from this observation that these matrix elements should be proportional to each 
other and calculations confirm this expectation. So, 

n s^s+.- = (^d\Y+Y? b \Q) = -y/2{dd\YP h tP b \ti) = -V2U ( i\d,u,b) . (13) 

Making the replacement u -B- d in Eq. (13), we obtain 

n s'^E-.+ = ^ u | E -go| ^ = y^2 (uu 0) = V2U?(u,d,b) . (14) 

In estimating the coupling constants of SSP, SAP and AAP, it is enough to consider the 
Tj b — > S°P, S fe ° — > E b P and S° — > E b P transitions, respectively. All remaining transitions 
can be obtained from these transitions with the help of the appropriate transformations 
among quark fields. Relations among the invariant functions of the charmed baryons can 
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easily be obtained by making the the replacement b — > c and adding to charge of each 
baryon a positive unit charge. 

Performing similar calculations, one can obtain rest of the required expressions from the 
correlation functions in terms of the invariant function III, involving ir, K and i] mesons 
describing sextet-sextet, sextet-antitriplet and antitriplet-antitriplet transitions. In the 
present work, we neglect the mixing between 77 and r] mesons. It should also be noted here 
that all coupling constants for the SSP, SAP and AAP are described only by one invariant 
function in each class of transitions, but the forms of the invariant functions in each group 
of transitions are different. 

The invariant function ITi responsible for the ->■ ££P, ->■ E^P and ~£ ->■ E° b P 
transitions can be calculated in deep Euclidean region, — p 2 — > +00 and — (p + q) 2 — > +00 
using the operator product expansion (OPE) in terms of the distribution amplitudes (DA's) 
of the pseudoscalar mesons and light and heavy quark operators. Up to twist-4 accuracy, 
the matrix elements (P(q) \q(x)Tq(0) \ 0) and (P(q) \q(x)G ^(0)1 0), where T is any arbitrary 
Dirac matrix, are determined in terms of the DA's of the pseudoscalar mesons, and their 
explicit expressions are given in [18-20]. 

The light and heavy quark propagators are calculated in [21], and [22], respectively. 
Using expressions of these propagators and definitions of DA's for the pseudoscalar mesons, 
the correlation function can be calculated from the QCD side, straightforwardly. Equating 
the coefficients of the structure ^75 of the representation of the correlation function from 
hadronic and theoretical sides, and applying the Borel transformation with respect to the 
variables p 2 and (p + q) 2 in order to suppress the contributions of the higher states and 
continuum, we obtain the following sum rules for the strong coupling constants of the 
pseudoscalar mesons with sextet and antitriplet baryons: 

(i)2 (i)2 

ff « = -^-^e^f + ^f , (15) 

A^A 2 

where, i — 1, 2 and 3 for sextet-sextet, sextet-antitriplet and antitriplet-antitriplet, re- 
spectively and M 2 and Mf are the Borel masses corresponding to the initial and the final 
baryons. Since the masses of the initial and final baryons are practically equal to each 
other, we take M\ = Mf = 2M 2 ; and A? and A^ } are the residues of the initial and final 
baryons, respectively, which are calculated in [23]. The explicit expressions for are 
quite lengthy and we do not present all of them here. As an example, we only present the 
explicit expression of the n^, which is given as: 



e m y M2 - m r/ M2 n?(u,d,b) = 



,,^ ;2 M'm%j v [l 2 - m 2 Q h]Uuo) + [ -^^M 4 m 2 Q ^ v { (i 3 (T, 1) - 2i 3 (T,vfjl 



- 2m 2 n 



(1-/3) 



128V27T 2 



i 3 (T, 1) - 2t 3 (T,v) + (1 - -ft, )^(txo)]/ 3 } 

M 2 m 2 v m Q f P \m 2 Q K(u )I 2 - 2(i 2 (V h 1) - 2i 2 (V ± , l^h 



' 2m 2 Q 



i 2 {A h 1) - 2(i 2 (V| h 1) - i 2 (V ± , 1) + i 2 (^|,- 
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?2\ 



' J -M 2 



3m^m 2 Q ii v {2rn 2 Q h - I^j (i 2 (T, 1) - 2i 2 (7») - <±(dd) fan 2 ^(uq) 



+ 



96 V / 27T ; 
(1 - (5 2 ) 

3M ^ M6 m r m Q m2 ofv(dd)A(u ) 

' j(<ta){(l - P 2 )m 2 v f P [5A(«o) + 12(W,|, 1) + i 2 (V| h 1) - 2i 2 (^||,u)) 



r4 " m Q m 0< 



2304 V^M 4 

- 8m Qf Pp(l - j%)[3 + 0(2 + 3/3)]0 CT ( Mo ) 

- 864 ^ m2 ^q(^){9(1 - P 2 )m Q f v (m 2 v A(u ) + mjfa„(uo)) + 2 [5 + 0(4 + 50)]m^(l - p^)0 CT K)} 

- {(1 - P 2 )fv{dd) [Qm 2 v A(u ) - 12m 2 P (i 2 (A h 1) + i 2 (V h 1) - 2i 2 {A h v)) + mg^(«o) 

+ 8[3 + /3(2 + 3/3)]m Q /i P (l - ^)(Jrf)0 a ( Mo )} (16) 
where J n is defined as: 



J ™ = / 2 rfs ^ 



and other parameters and functions as well as the way of continuum subtraction are given 
in [14]. To shorten the equation, we have ignored the light quark masses as well as terms 
containing gluon condensates in the above equation, but we take into account their contri- 
butions in numerical analysis. 

3 Numerical results 

In this section, we present the numerical results of the sum rules for strong coupling con- 
stants of pseudoscalar mesons with sextet and antitriplet heavy baryons, which are obtained 
in the previous section. The main input parameters of LCSR are DA's for the pseudoscalar 
mesons which are given in [18-21]. The other input parameters entering to the sum rules are 
(qq) = -(0.24±0.001) 3 GeV 3 , m 2 = (0.8±0.2) GeV 2 [24], / w = 0.131 GeV, f K = 0.16 GeV 
and = 0.13 GeV [18]. 

The sum rules for the SSP, SAP and AAP coupling constants have three auxiliary 
parameters: Borel mass parameter M 2 , continuum threshold Sq and the arbitrary parameter 
[5 which exists in the expression in the expression for the interpolating currents. Obviously, 
the result for any measurable physical quantity, being coupling constant in the present case, 
should be independent on them. Therefore, our primary goal is to find such regions of these 
parameters, where coupling constants exhibits no dependence. 

The upper limit of M 2 is determined by requiring that the continuum and higher states 
contributions should be small compared to the total dispersion integral. The lower limit 
can be obtained from the condition that the condensate terms with highest dimensions 
contributes smaller compared to the sum of all terms. These two conditions leads to the 
working region, 15 GeV 2 < M 2 < 30 GeV 2 for the bottom baryons and 4 GeV 2 < M 2 < 
12 GeV 2 for the charmed ones. The continuum threshold is not totally arbitrary but it 
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depends on the energy of the first excited state with the same quantum numbers as the 
interpolating current. We choose it in the domain between s = (mg + 0.5) 2 GeV 2 and 
so = {jriB + l) 2 GeV 2 . As an example, let us consider the S 6 ° — y S 6 °7r° transition. In Fig. 
(2), the dependence of the strong coupling constant for the E'® — y S^ 7r° transition on M 2 is 
considered at different fixed values of f3 and a fixed value of s . We observe from this figure 
that the coupling constant has a good stability in the "working region" of M 2 . In Fig. (3), 
we present the dependence of the strong coupling constant for the S b ° — y S 6 °7r° transition on 
cos 6* at several fixed values of s and at M 2 = 22.5 GeV 2 , where the angle 9 is determined 
from j3 = tan#. From this figure, we see that the dependence of the coupling constant 
on s diminishes when the higher values of the continuum threshold are chosen from the 
considered working region. From this figure, we also observe that the strong coupling 
constant for the —y S fe °7r° decay becomes very large near the end points (cos# = ±1) 
and have zeros at some finite values of the cos 9. This behavior can be explained as follows. 
From Eq. (15) we see that the coupling constant is proportional to (i * w 11^. In general, 

zero's of the nominator and denominator does not coincide since the OPE is truncated. 
In other words, calculations are not exact. For this reason, these points and any region 
between them are not reliable regions for determination of physical quantities and suitable 
regions for cos# should be far from these regions. It follows from Fig. (3) that in the 
region, —0.5 < cos 9 < +0.3, the coupling constant seems to be insensitive to the variation 
of cos 6*. Here, we should also stress that our numerical results lead to the working region, 
—0.6 < cos# < +0.5 common for masses of all heavy spin 1/2 baryons, which includes the 
working region of cos 9 for the coupling constant. This region lie also inside the more wide 
interval of cos 9 obtained from analysis of the masses of the non strange heavy baryons in 
[10,25,26]. In general, the working region of cos 9 for masses and coupling constants can 
be different, but in some cases as occur in our problem these regions coincide. 

Similar analysis for the strong coupling constants of the light pseudoscalar mesons with 
sextet and antitriplet heavy baryons are performed and the results are presented in Tables 
(2), (3) and (4). In these Tables we also present the predictions for the coupling constants 
coming from the Ioffe currents when (5 = — 1. The errors in the values of the coupling 
constants presented in the Tables (2), (3) and (4) include uncertainties coming from the 
variations of the s , f3 and M 2 as well as those coming from the other input parameters. 
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Table 2: The values of the strong coupling constants g for the transitions among the sextet 
and sextet heavy baryons with pseudoscalar mesons. 
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Table 3: The values of the strong coupling constants g for the transitions among the sextet 
and antitriplet heavy baryons with pseudoscalar mesons. 

We see from these Tables that, there is substantial difference between the predictions of 
the general current and the Ioffe current, especially for the strong coupling constants of the 
antitriplet-antitriplet heavy baryons with pseoduscalar mesons, which can be explained as 
follows. As a result of the analysis of the dependence of the coupling constants on cos 9 we 
see that the value /3 — — 1 belongs to the unstable region. Therefore, a prediction at this 
point of P is not reliable. 

Finally, we compare our results with those existing in literature. In various works, the 
coupling constant S c — > A c n is estimated to be 

{8.88, [27] (relativistic three-quark model) , 
6.82, [28] (light-front quark model), 
10.8 ±2.2, [9] (LCSR) , 
6.5 ± 2.4, (our result) (LCSR) . 
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Table 4: The values of the strong coupling constants g for the transitions among the 
antitriplet and antitriplet heavy baryons with pseudoscalar mesons. 



We see that, within errors our result is close to the results of [9,27,28]. The coupling 
constant for the EqEqtt transition LCSR is estimated to have the values g-^-^ = 1.0±0.5 
and g =f >~^ b7r = 1.6 ±0.4, which are slightly larger compared to our predictions. Finally, the 
coupling constant g s c->s C 7r j g ca J cu l a ted in [9] and it is obtained that g E c->s C 7r = — 8.0 ± 1.7, 
which is in quite a good agreement with our prediction. 

In conclusion, the strong coupling constants of light pseudoscalar mesons with sextet and 
antitriplet heavy baryons are studied within LCSR. It is shown that, all coupling constants 
for the sextet-sextet, sextet-antitriplet and antitriplet-antitriplet transitions are described 
by only one invariant function in each class. 
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Appendix A : 



Here in this appendix, we present the expressions of the correlation functions in terms 
of invariant function II ^ involving tt, K and r] mesons. 



• Correlation functions describing pseudoscalar mesons with sextet-sextet baryons. 
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Correlation functions responsible for the transitions of the sextet-antitriplet baryons. 
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1 

^6 L 
1 
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Uf\u,s,b) + 2nj 2) (s,M,6) 

nS 2) (rf,s,6) + 2ni 2) (s,rf,6) 
nS 2) (M,rf,6) -nS 2) (rf,w,6) 



Correlation functions appearing in the antitriplet-antitriplet pseudoscalar meson tran- 
sitions. 



V2IF- 



nS 3) ( M ,s,6) 
nj 3) (d,s,&) 



nS 3) ( M ,rf, 6) -n^ j (rf,M,6) 



(3), 



v/2 

ni 3) («,«,6) , 

-Tl?(u,d,b) , 



10 



j-jH-^-S-j?! = J_ 



ni 3) (M, s ,6) -2ni 3) (s,w,6) 

II< 3) (d,s,&) -2U ( l\s,d,b) 



Uf\d,u,b) + Uf\u,d,b) 



The expressions for the charmed baryons can easily be obtained by making the replace- 
ment b — )• c and adding to charge of each baryon a positive unit charge. 
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Figure 2: The dependence of the strong coupling constant for the S^ — y E'®ir° transition 
at several different fixed values of /3, and at s = 40.0 GeV 2 . 




Figure 3: The dependence of the strong coupling constant for the — y ^tt° transition 
on cos 9 at several different fixed values of s , and at M 2 = 22.5 GeV 2 . 
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